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FOREWORD 

This  paper  was  presented  by  Mr.  W.  G.  Shockley  at  the  Fifth  Meeting 
of  the  Quadripartite  Standing  Working  Group  on  Ground  Mobility,  which 
was  held  in  Kingston  and  Parry  Sound,  Ontario,  Canada,  during  9-19  August 
1965.  The  information  contained  herein  summarizes  studies  conducted  for 
the  U.  S.  Array  Engineer  Waterways  Experiment  Station  (WES)  by  Chrysler 
Corporation  under  Contract  No.  DA-22- 079-eng-403  and  by  FMC  Corporation 
under  Contract  No.  DA-22-079-eng-4ll. 

The  Chrysler  study  was  a  portion  of  the  Mobility  Environmental  Research 
Study  (MERS),  sponsored  by  the  Office,  Secretary  of  Defense,  Advanced  Re¬ 
search  Projects  Agency,  Directorate  of  Remote  Area  Conflict,  under  Order 
No.  400.  WES  is  the  prime  contractor  for  the  MERS  project  and  the  U.  S. 
Army  Materiel  Command  ( AMC )  is  the  service  agent.  The  Chrysler  study  was 
reported  in  "A  Study  of  the  Vehicle  Ride  Dynamics  Aspects  of  Ground 
Mobility, WES  Contract  Report  No.  3“ll4,  volumes  1  and  2  dated  March  1965 
and  volumes  3  and  4  dated  April  1965. 

The  FMC  study  was  conducted  in  furtherance  of  Department  of  the  Army 
Research  and  Development  Projects  1-V-0-25001-A-131,  Military  Evaluation 
of  Geographic  Areas,  which  is  sponsored  by  AMC.  The  study  was  reported 
in  "A  Computer  Analysis  of  Vehicle  Dynamics  While  Traversing  Hard  Surface 
Terrain  Profiles,"  WEC  Contract  Report  No.  dated  February  1966. 

Colonel  John  R.  Oswalt,  Jr.,  CE,  was  Director  of  the  WES  at  the  time 
of  the  preparation  and  presentation  of  this  paper.  Mr.  J.  B.  Tiffany 
was  Technical  Director.  . 


54048 


BUMPS  AND  GRINDS:  STUDIES  IN  BODY  MOTION 
by 

W.  G.  Shockley* 

Introduction 

1.  About  four  years  ago,  when  the  Waterways  Experiment  Station  (WES) 
made  its  first  tentative  steps  toward  the  development  of  a  quantitative 
cross-country  mobility  prediction  system  for  vehicles,  a  series  of  cross¬ 
country  speed  tests  was  conducted  at  the  Yuma  Test  Station  in  the  American 
desert.  Among  other  things,  these  tests  revealed  that  small-scale  surface 
geometry  features  were  among  the  most  important,  if  not  the  most  important, 
of  the  terrain  attributes  controlling  the  cross-country  speed  of  military 
vehicles.  It  was  found  that  the  vehicles  rarely  attained  the  speed  permitted 

by  the  strength  of  the  surface  they  were  crossing;  instead,  they  were  continually 
either  accelerating  away  from  a  surface  geometry  feature  or  decelerating  in 
preparation  for  encountering  one.  It  was  obvious  from  the  tests  that  one  of 
the  most  critical  vehicle  characteristics  was  "ride"  quality;  that  is,  the 
dynamic  response  of  the  machine  to  a  vertical  obstacle  (i.e.  a  surface  feature 
which  results  in  vehicle  motion  in  the  vertical  plane  as  the  vehicle  surmounts 
it). 

2.  As  a  result  of  this  experience,  the  WES  began  an  investigation  of 
the  state  of  the  art  with  respect  to  ride  dynamics,  since  it  was  quite  obvious 
that  any  satisfactory  mat  .:r’i atical  model  of  terrain-vehicle  relations  would 
have  to  include  the  suspension  characteristics  of  the  vehicle.  Our  investiga¬ 
tions  revealed  that  the  automotive  industry  and  such  laboratories  as  the  Land 
Locomotion  Laboratory  had  been  working  on  the  general  problem  for  some  time, 
and  that  they  had  made  great  strides  in  both  the  measurement  gold  evaluation 
of  ride  quality.  Hie  approach  taken  by  these  organizations  was  that  "ride" 
quality  was  a  function  of  a  spectrum  of  vibrations  (i.e.  small  amplitude  but 
relatively  high  frequency  cyclic  motion),  and  therefore  the  analytical  approach 
was  to  treat  both  terrain  input  and  dynamic  vehicle  response  in  terms  of  wave 
motion.  Hie  most  useful  device  was  to  reduce  both  terrain  characteristics  and 
vehicle  response  to  power  spectral  density  (PSD)  relations,  and  then  correlate 
the  two  statements  by  involving  a  jury. 

*  Engineer,  Chief,  Mobility  and  Environmental  Division,  U.  S.  Army  Engineer 
'Waterways  Experiment  Station,  Vicksburg,  Miss. 


3.  When  we  examined  this  problem  in  the  light  of  our  experiences 
in  cross-country  mobility,  a  school  of  thought  arose  which  questioned 
whether  this  approach  could  provide  an  adequate  evaluation  of  cross¬ 
country  vehicles.  Experience  seemed  to  indicate  that  speed  was  limited 
by  the  magnitudes  of  the  bumps  which  were  generated  as  the  vehicle* 
surmounted  discrete  vertical  obstacles,  rather  than  by  the  cumulative 
grind  of  vibration.  One  of  the  mathematicians  working  on  related  matters 
(Dr.  James  Dugundji,  University  of  Southern  California)  insisted  that  no 
form  of  vibrational  analysis  would  elucidate  the  cause-effect  relations 
vhen  a  vehicle  hits  a  bump.  He  found  a  substantial  body  of  support  at 
the  WES,  and  as  a  result  we  are  currently  favored  by  at  least  three 
schools  of  thought  on  this  question. 

k.  The  first  school  is  convinced  that  vibrational  analysis  is  an 
adequate  approach,  because  it  is  held  that  the  response  of  a  vehicle  to 
a  discrete  vertical  obstacle  is  simply  a  high-frequency,  low-amplitude 
vibration,  and  is  therefore  susceptible  of  treatment  by  PSD  analysis. 

The  second  school  insists  that  design  improvement  in  suspension  systems 
is  dependent  upon  the  ability  to  follow  the  force  systems  through  the 
individual  vehicle  components,  and  that  therefore  one  needs  a  form  of 
analysis  based  on  fundamental  force  systems.  That  is,  one  must  be  able 
to  predict  what  a  unique  surface  feature  will  produce  in  the  way  of 
vehicle  response.  The  third  school  straddles  the  fence;  it  argues  that 
the  PSD  analysis  is  entirely  adequate  for  small-amplitude  surface  irregu¬ 
larities,  but  that  the  discrete  prediction  scheme  will  have  to  be  used 
for  large  obstacles, 

5.  The  question  between  bumps  and  grinds  coula  not  be  resolved  by 
debate,  so  it  was  decided  to  investigate  the  potentiality  of  both.  The 
internal  capabilities  at  the  WES  were  inadequate  to  the  task,  so  the 
problem  was  turned  over  to  two  contractors:  the  FMC  Corporation  is 
investigating  the  "bump"  approach;  the  Chrysler  Corporation  is  extending 
its  previous  work  on  the  "grind”  approach. 

Runps:  The  SMC  Corporation  Study  of  Ride  Dynamics 

6.  The  point  of  departure  for  this  study  was  the  desire  to  develop 
a  completely  objective  method  of  predicting  the  absolute  displacements 
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generated  at  any  position  in  a  vehicle  as  a  response  to  an  encounter 
with  any  surface  geometry  feature  at  any  speed.  The  implications  of 
this  are: 

a.  That  an  accurate  mathematical  model  of  the  surface  be 
esployed  as  terrain  input. 

b.  That  the  mathematical  model  of  the  vehicle  be  a  component- 
by-component  construction,  so  that  information  on  displacement  can  be 
recovered  at  any  point. 

7.  Any  attempt  to  mathematically  model  a  terrain  for  the  purpose 
of  predicting  vehicle  dynamic  response  raises  two  questions: 

a.  Can  the  model  be  made  to  fit  the  real  terrain  within 
acceptable  limits  of  accuracy? 

b.  Will  the  analytically  predicted  response  of  the  vehicle 
traversing  this  terrain  agree  within  acceptable  limits  of  accuracy  to 
reality? 

8.  Because  there  existed  no  prior  data  on  the  degree  of  fit  required 
by  item  a  in  paragraph  7,  it  was  decided  to  model  the  aurface  geometry 

by  fitting  a  Fourier  series  to  it.  The  advantage  is  simply  that  the 
accuracy  of  fit  is  a  function  of  the  number  of  terms  in  the  Fourier 
series;  greater  or  lesser  degrees  of  fit  can  be  achieved  by  adding  or 
subtracting  terms  from  the  Fourier  expression.  Therefore,  a  Fourier 
series  was  used  having  the  general  form: 


i*l 

where:  X  =  microgeometry  evaluation  at  horizontal  displacement  y 

&0*  twice  the  average  value  of  X 
bi=  sine  coefficient  for  the  i  term 
phase  angle  for  the  i  term 
i  ■  series  term  number 
p  =  fundamental  period  length 
n  *  total  number  of  terms  in  series 
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This  sine-phase  angle  form  was  preferred  over  the  sine-cosine  form 
because  computer  time  was  saved  in  evaluating  the  trigonometric  func¬ 
tions.  A  pure  sine  or  cosine  series  could  not  be  used  since  both  odd 
and  even  functions  were  components  of  the  overall  curve  fit.  This 
property  was  introduced  because  a  method  for  extending  the  profile 
had  to  be  employed}  the  Fourier  series  describes  a  cyclic  function 
extending  from  +  «•  to  -  •,  while  the  surface  geometry  profile  is  of 
finite  length.  The  method  chosen  was  to  repeat  the  given  profile 
indefinitely,  so  that: 

f(y  +  np)  =  f(y)  (2) 

as  illustrated  in  fig.  1. 

9.  It  was  also  necessary  to  make  the  start  and  end  of  the  profile 
at  the  same  level,  f(np)  =  f(0),  so  that  an  abrupt  change  would  not 
occur  at  the  end  of  one  cycle  and  the  start  of  another.  This  level  was 
selected  at  zero,  f(np)  *  0,  so  that  the  initial  conditions  of  a  vehicle 
traversing  the  terrain  could  all  be  set  equal  to  zero. 

10.  The  digital  computer  program  for  evaluating  the  Fourier  series 
terms  to  fit  this  microrelief  is  presented  in  Appendix  A.  This  program 
was  designed  to  accept  surface  geometry  profile  data  as  they  come  in 
from  the  field.  It  had  been  found  that  the  surveying  method  which  was 
most  conservative  of  time  and  effort  was  to  locate,  in  terms  of  distance 
and  elevation,  each  "break"  in  slope.  On  curved  surfaces,  many  such 
points  would  be  selected;  on  planar  surfaces,  the  points  would  be  far 
apart.  Thus,  the  surveying  results  consist  essentially  of  a  series  of 
irregularly  spaced  x-y  coordinates.  The  existing  program  accepts  these 
sets  of  coordinates  locating  the  points  of  discontinuity  and  calculates 
linear  functions  between  them.  These  linear  functions  are  used  as  the 
f(y)  for  the  Fourier  integrations.  Up  to  100  coordinate  sets  may  be 
specified.  The  program  output  consists  of  the  constant  Sq,  the  sine 
coefficients  bi,  and  the  phase  angles  a^  where  i  is  the  series  term 
number.  The  total  number  of  terms  may  be  any  number  desired  up  to  52. 
Calculated  amplitudes  are  compared  with  the  input  coordinates  and  the 
error  of  fit  is  included  in  the  output. 
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11*  As  a  test  of  the  program,  a  simple  s.rface  geometry  feature 
(fig.  2)  was  used  to  check  the  program.  A  Uo-term  Fourier  series  fit 
the  bump  profile  with  a  maximum  error  of  1.5*  (i.e.  the  maximum  depar¬ 
ture  of  the  computed  profile  from  the  real  profile  was  1.5*  of  the  total 
height  of  the  real  obstacle)*  in  addition,  a  number  of  surface  geometry 
configurations  extracted  from  Appendix  E  of  WES  Technical  Report  5-625  # 
were  fitted  with  a  52-term  Fourier  series.  The  maximum  fit  error  for  each 


microgeometry 

feature  was: 

Percent 

Illustrated  on 
fig.  No. 

Site 

!52: 

Rice-field  bund,  A 

1.68 

3 

Rice- field  bund,  B 

1,1*6 

Field  prepared  for 
unidentified  row  crop 

2.66 

1* 

Field  surface  fallow 
after  ride 

6.6h 

Site 

1 93: 

Rice- field  bund 

2.1*5 

Site 

2155 

Canalized  water  channel 

8.19 

5 

Site 

U5i: 

Stream  channel 

6.02 

6 

Site 

532: 

Boulder-strewn  slope,  A 

fc.57 

Boulder-strewn  slope,  B 

2.07 

7 

Site 

92i,: 

Terraced  rice  field,  A 

2.25 

Terraced  rice  field,  B 

1». 70 

These  profiles  were  plotted  by  scaling  the  plates  in  TR  5-625.  The 
terrain  heights  listed  under  "input"  are  the  scaled  values;  the  "output” 
heights  are  the  corresponding  values  of  the  Fourier-series  curve.  A 
comparison  of  these  two  height  columns  indicate  the  absolute  fit  achieved. 
The  accuracy  of  fit  was  considered  satisfactory  for  all  profiles. 

12.  The  next  step  in  the  process  was  to  develop  a  three-dimensional 
mathematical  model  of  a  vehicle  (an  M37  3A-ton  truck  was  used  as  an 
example)  and  test  it  using  the  Fourier  expression  as  the  surface  geometry 

U.  S.  Army  Engineer  Waterways  Experiment  Station,  CE.  Environmental 

ghgtors  Affecting  Ground  Mobility  in  Thailand.  Technical  Report - 

Ao.  5-625  with  Appendices  A-H  (Vicksburg,  Mississippi,  May  1963) 


input.  The  truck  body  was  allowed  three  degrees  of  freedom  (bounce, 
pitch,  an  ;,  roll)  and  was  supported  by  four  spring-damper  combinations 
■  represe.-  .~g  the  suspension.  The  suspension  was  supported  by  a  differ¬ 
ential  assembly  front  and  rear,  each  having  t.  v  degrees  of  freedom 
(bounce  and  roll).  The  differential  assemblies  were  supported  by 
spring-daape  combinations  representi.i^  the  tires.  This  model  with  the 
necessary  symbols  is  shown  in  fig.  8.  Equations  for  input  parameters 
(tire  deflections  and  deflection  rates,  and  suspension  deflections  and 
4: 'I? action  rates)  are  given  in  table  1.  The  model  requires  seven  simul- 
iane.s  .  cond-order  ordinary  differential  equations  with  first  deriva- 
;i’  /*..  cse  equation®  are  given  in  table  2. 

* ’ .  A  preliminary  computer  program  was  developed  to  solve  the  seven 
differential  equations  of  motion.  This  initial  program  placed  no  restraint 
on  suspension  deflection  and  assumed  that  the  tires  remain  in  contact  with 
the  ground.  These  conditions  are  probably  realistic  for  slow  travel  and 
smal...  responses.  A  Fourier  series  description  of  the  surface  geometry  is 
used  for  the  displacement  forcing  function  with  the  option  of  different 
time  lags  for  each  wheel  to  simulate  the  traversing  of  the  terrain  at  an 
angle  to  the  ridge  lines,  thereby  inducing  a  roll  motion  i.  the  response. 
The  program  output  consists  of  the  displacement,  velocity,  and  accelera¬ 
tion  of  the  truck  body  in  bounce,  pitch,  and  roll.  Figs.  9  through  lU 
illustrate  the  predictions  of  bounce  and  pitch  displacements,  velocities, 
and  accelerations.  Hie  response  curves  derived  from  th~  standard  two- 
dimensional,  8 ingle-mass  mathematical  model  which  had  previously  been 
used  are  included  for  comparison. 

lfc.  The  next  logical  step  consists  of  two  parts: 

a.  Refine  the  three-dimensional  model  to  take  care  of  such 
things  as  wheels  leaving  the  ground,  the  suspension  system  "bottoming 
out,"  and  the  nonlinear  deflections  and  recoveries  of  suspension  system 
elements;  and 

b.  Test  the  model  against  real  performance  values  obtained  in 
the  field. 

15.  Both  of  these  phases  have  been  completed.  The  model  is  now 
about  as  elegant  as  it  is  practical  to  make  it  at  the  present  time;  it 


already  needs  a  computer  of  rather  phenomenal  capacity.  A  diagrammatic 
mathematical  model  of  a  generalized  wheeled  vehicle  (only  very  slight 
modifications  are  required  to  accommodate  tracks)  is  illustrated  in 
fig.  15.  Furthermore,  the  model  was  altered  to  include  such  things  as: 

Variable  rolling  radii,  as  illustrated  in  fig.  l6 

l 

b.  Horizontal  ami  vertical  force  reactions  at  the  wheel  ground 
contact  point**  (fig.  16) 

c.  Wheel  bounce,  so  that  it  loses  contact  with  the  surface 

(fig.  17) 

d.  Nonlinear  and/or  discontinuous  suspension  spring  constants, 
to  take  into  account  the  situations  in  which  the  vehicle  "bottoms  out" 
(fig.  17) 

e.  Nonlinear  and  discontinuous  responses  of  shock  absorbers 

(fig.  18) 

f.  Variable  horizontal  speed. 

The  new  equations  describing  body  bounce,  body  pitch,  and  body  roll  are 
presented  in  table  3.  This  new  model  was  then  "run"  over  two  of  the 
terrain  types  previously  used  for  evaluating  the  simplified  model  (figs.  1* 
and  7)  and  the  predicted  motions  (i.e.  roll,  pitch,  bounce,  and  horizontal 
speed)  recorded. 

16.  Next,  models  of  these  terrain  features  were  constructed 

of  reinforced  concrete,  in  two  different  versions.  The  first  was  so 
arranged  that  the  obstacles  were  encountered  at  right  angles  to  the  axis 
of  motion  (fig.  19),  and  the  second  version  was  so  arranged  that  the 
angle  of  impact  was  60°  (fig.  20). 

17.  The  test  vehicle,  an  M37»3A-ton  UxU  truck,  was  disassembled 
and  all  of  the  required  parameters  were  measured  in  the  laboratory,  to 
make  certain  that  the  vehicle  descriptors  were  correct.  The  vehicle  was 
then  put  back  together,  appropriately  instrumented,  and  run  through  the 
courses  under  extremely  carefully  controlled  conditions,  with  the  motions 
being  recorded  continuously  by  two  different  recording  schemes.  One 
employed  vehicle-mounted  accelerometers ;  the  other  used  continuous  high¬ 
speed  photography. 
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18.  The  results  are  still  being  analyzed,  but  it  is  obvious  even 
at  the  present  time  that  the  predicted  responses  and  measured  responses 
fit  with  remarkable  fidelity.  The  maximum  error  is  about  8J,  well  within 
the  1C%  error  which  had  been  accepted  as  tolerable.  Fig.  21  illustrates 
a  portion  of  one  of  the  response  curves. 

19.  If  it  be  accepted  that  the  mathematical  model  truly  predicts 
dynamic  vehicle  responses  within  an  acceptable  level  of  accuracy,  then 
it  is  clear  that  a  powerful  tool  for  assisting  in  the  design  of  cross¬ 
country  vehicles  is  in  hand.  It  is  anticipated,  for  example,  that  a 
designer  will  be  able  to  vary  response  characteristics  or  other  vehicle 
parameters  one  at  a  time,  and  evaluate  the  results  on  the  machine  as  a 
whole.  It  seems  obvious  that  design  solutions  very  close  to  optimum 
could  be  achieved  relatively  quickly  and  cheaply  by  this  method. 

Grinds:  The  Chrysler  Corporation  Investigations  of  RLde  Dynamics 
Introduction 

20.  Up  to  this  point  the  WES  activities  concerned  with  the  effects 
of  singular  vertical  obstacles  which  cause  a  transient  response  of  the 
vehicle  system  have  been  described.  In  this  case,  vehicle  speed  is 
primarily  limited  by  the  operator's  anticipation  of  the  "bump”  which  will 
result  when  his  vehicle  strikes  the  obstacle. 

21.  In  addition  to  these  studies,  WES  has  conducted  as  part  of  the 
Mobility  Environmental  Research  Study  (MERS)  being  conducted  for 
Advanced  Research  Projects  Agency  (ARPA)  a  study  which  is  intended  to 
elucidate  the  relations  between  terrain  roughness  and  vehicle  vibration. 
This  study  has  been  accomplished  under  contract  by  Defense  Engineering, 
Chrysler  Corporation,  with  Mr.  B.  D.  Van  Ceusen  as  the  principal  investi¬ 
gator.  The  criterion  for  vehicle  evaluation  is  that  of  human  response  to 
vehicle  vibration;  hence  the  "grinds.”  In  this  case,  it  is  assumed  that 
the  vehicle  operator  will  vary  speed  until  the  vibration  of  the  vehicle 
(i.e.  the  "ride")  is  tolerable  or  acceptable.  Factors  which  may  determine 
sustainable  speed  include  the  ability  of  the  operator  to  control  the 
vehicle,  human  tolerance  to  vibration,  and  damage  to  or  destruction  of 
the  vehicle  or  cargo. 
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22,  The  purpose  of  the  Chrysler  Corporation  vehicle  ride  dynamics 
study  was  to  assess  the  adequacy  of  present  analytical  systems  for  pre¬ 
dicting  the  ride  performance  of  military  vehicles  operating  cross-country 
over  irregular  surface  profiles  at  speeds  up  to  25  mph,  and  to  establish 
areas  that  need  further  development  to  effect  good  correlation  between 
the  analytical  prediction  and  actual  vehicle  behavior.  The  study  was 
basically  a  state-of-the-art  determination  of  analytical  systems  appli¬ 
cable  to  the  problem  of  vehicle  ride  dynamics.  As  stated  earlier,  the 
dtudy  was  not  concerned  with  shock  or  impact  aspects  of  vehicle  dynamics 
caused  by  singular  vertical  obstacles  but  with  the  more  or  less  steady- 

4 

state  vehicle  vibration,  the  subjective  quality  of  which  is  termed  "ride" 
in  the  literature.  A  limited  amount  of  effort  was  expended  in  field 
experimentation  to  develop  new  or  improve  existing  analytical  systems 
from  state-of-the-art  knowledge. 

23.  A  review  of  pertinent  literature  indicated  that  some  qualitative 
criteria  have  been  developed  which  will  allow  approximate  measurements 

of  vehicle  ride.  However,  it  was  recognized  that  a  definite  need  exists 
for  quantitative  analytical  procedures  which  would  permit  designers  to 
predict  ride  performance  of  a  vehicle  over  a  given  irregular  surface,  as 
a  function  of  vehicle  speed.  A  minimum  of  three  elements  of  the  problem 
must  be  attacked  in  order  to  obtain  numerical  descriptors  as  ridi  indices 
for  vehicle  comparison.  First,  it  is  necessary  to  quantify  t.he  subjective 
reactions  of  humans  in  order  to  obtain  a  meaningful  scale  describing  ride 
quality.  Second,  mathematical  models  of  terrain-vehicle  interaction  must 
be  developed,  so  that  the  vibration  response  within  the  vehicle  can  be 
associated  with  a  particular  terrain  input.  Third,  a  practical  method 
of  obtaining  and  analyzing  terrain  profiles  must  be  developed.  When 
these  elements  of  the  problem  have  been  resolved,  it  will  be  possible  to 
make  design  decisions  from  computer  model  studies,  since  it  will  be 
possible  to  evaluate  the  effect  of  tradeoffs  between  ride  and  other 
vehicle  design  considerations  such  as  cost,  weight,  complexity,  and 
reliability. 
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Human  response  to  vehicle  vibration 

2U,  One  criterion  of  interest  for  vehicle  evaluation  is  that  of 
jnuman  response  to  vibration.  Thus,  criteria  for  defining  human  response 
are  necessary  in  order  to  determine  the  nature  of  the  output  of  terrain- 
vehicle  models.  It  is  obvious  that  the  form  of  the  output  will  to  some 
extent  determine  the  form  of  the  model  which  is  used.  This,  in  turn,  will 
determine  to  some  extent  the  form  of  the  terrain  measurements  used 
as  inputs  to  the  model. 

25.  A  number  of  previous  research  programs  have  attempted  to  define 
human  response  to  vibration.  In  most  of  these  programs,  ride  simulators 
with  sinusoidal  inputs  have  been  used  and  attempts  made  to  correlate 
comfort  with  frequency.  While  this  information  does  not  directly  pertain 
to  the  type  of  vibration  associated  with  most  cross-country  vehicles,  a 
good  deal  of  background  can  be  obtained  from  study  of  the  literature  de¬ 
scribing  these  experiments.  Three  problems  exist  in  the  interpretation 
of  experiments  with  ride  simulators  using  sinusoidal  inputs  in  the  light 
of  the  present  objective: 

a.  The  determination  of  the  relationship  (if  any)  between 
sinusoidal  inputs  and  random  vibration,  such  as  that  excited  by  terrain 
irregularities.  It  is  suspected  that  a  sine  wave  of  a  single  frequency 
is  more  objectionable  to  the  human  than  a  random  input  having  the  same 
vibratory  energy. 

b.  Most  shake  table  experiments  have  involved  vibration  in  a 
single  plane.  It  is  strongly  suspected  that  vibration  imposed  simul¬ 
taneously  in  more  than  one  mode  will  not  be,equivalf nt  to  the  sum  of 
responses  to  each  mode  individually. 

c.  The  criteria  used  for  determining  equacorafort  levels  in 
most  shake  table  experiments  are  vague  and  subjective.  The  usual  de¬ 
scriptive  terms  are  mildly  annoying,  alarming,  unpleasant,  intolerable, 
etc.  It  is  difficult  to  interpret  vehicle  vibrational  activity  in  terms 
of  these  subjective  criteria,  since  different  people  respond  in  sub¬ 
stantially  different  ways.  It  is  suggested  that  a  numerical  scale  of 
some  form  is  necessary  to  determine  meaningful  differences  between 
vibrational  levels. 
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26.  Fig.  22  shows  a  composite  of  eight  different  ride  simulator 
results  for  vibration  tolerance.  There  are  three  conclusions  that  can 
be  drawn  from  this  graph: 

a.  There  is  a  large  discrepancy  in  discomfort  level,  which  is 
attributed  to  the  vague  subjective  criteria  used.  Interpretation  of 
these  results  has  been  confounded  by  the  large  variety  of  experimental 
conditions  used  by  the  various  investigators. 

b.  A  plot  of  acceleration  versus  frequency  at  a  given  comfort 
level  is  relatively  flat.  It  seems  probable  that  the  derivative  of 
interest  is  the  acceleration  and  not  the  displacement,  velocity,  or  Jerk 
{time  rate  change  of  acceleration). 

c.  There  is  an  area  between  5  and  7  cps  where  the  human  is  more 
sensitive  to  vibration.  This  sensitivity  is  attributed  to  visceral 
resonance  or  resonance  of  the  abdominal  cavity.  Vibration  in  this  fre¬ 
quency  range  is  quite  disturbing,  and  the  visceral  resonance  range  is 
easily  distinguishable  in  a  ride  simulator  which  is  driven  slowly  through 
this  frequency  range. 

27.  In  order  to  determine  the  feasibility  of  measuring  magnitudes 
of  accelerations  in  such  a  way  that  they  could  be  correlated  with  sub¬ 
jective  Judgments  of  ride  quality,  five  people  vere  driven  over  21 
separate  paved  roads  and  unimproved  trails  at  the  Chrysler  Chelsea 
Proving  Grounds  in  two  different  wheeled  vehicles  (M37,  l*xU  3/^-ton 
truck,  and  the  XMl*10,  8x8  2-1/2-ton  truck).  Judgments  of  the  ride 
quality  were  obtained  from  each  person  for  each  road  section  in  each 
vehicle.  Time  recordings  of  the  vibratory  accelerations  along  three 
perpendicular  axes  were  obtained  at  the  seat  for  each  vehicle  for  each 
road  section.  These  acceleration  traces  were  recorded  on  magnetic  tape 
and  converted  with  an  analog  computer  into  statistical  quantities  suitable 
for  analysis.  The  subjective  Judgments  were  compared  with  the  physical 
quantities  derived  from  the  acceleration  recordings.  This  correlation 
constituted  the  major  objective  of  this  investigation  and  was  performed 
vith  a  digital  computer. 
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28.  Table  A  Is  a  summary  of  some  of  the  correlation  coefficients 
from  this  study. 


Table  A 


Correlation 

Coefficients 

Between 

Ride  and  RMS  Acceleration  1 

Criterion 

M37 

XMtilO 

Both 

Vehicles 

Vertical 

.683 

.830 

.77k 

Dieckmann 

Transverse 

.803 

.8U9 

.1*26 

Fore  &  Aft 

.579 

.887 

.619 

Vertical 

.1*83 

.895 

.787 

Van  Deusen 

Transverse 

.377 

.817 

.M7 

Fore  &  Aft 

.393 

.880 

,6»*o 

2p.  Table  A  shows  correlation  coefficients  between  ride  and  root 
mean  square  (rms)  acceleration.  Two  different  criteria  are  compared  in 
this  table.  The  Dieckmann  criterion  emphasized  low  frequencies  and  the 
Van  Deusen  criterion  emphasized  the  area  associated  with  the  visceral 
resonance.  The  largest  correlation  coefficient  found  in  this  study  was 
.895  for  the  XMUlO  using  the  Van  Deusen  criterion.  There  is  a  consider¬ 
able  difference  between  the  correlation  found  in  the  M3T  vehicle  and  the 
XMUlO  vehicle.  The  frequency  distribution  is  determined  by  the  vehicle 
configuration,  i.e.  the  road  or  terrain  input  has  no  predominant  frequency 
components  and,  therefore,  the  natural  resonance  frequencies  of  the 
vehicle  determine  the  frequency  content  of  the  vibration.  This  is  probably 
the  reason  for  the  rather  large  difference  in  correlation  coefficients  for 
the  two  vehicles.  It  is  suggested  that  a  single  number  (properly  filtered 
rms  acceleration)  derived  from  actual  acceleration  measurements  may  be  a 
meaningful  ride  index  for  a  given  vehicle,  but  to  compare  ride  for  a 
number  of  vehicle  configurations,  it  is  necessary  to  have  more  than  one 
number  to  characterize  the  vibration. 

30.  Conclusions  reached  from  the  review  of  the  state  of  the  art 
and  the  feasibility  study  are  as  follows: 

a.  It  is  clear  tl  at  the  problem  of  correlating  acceleration 
values  with  subjective  ride  sensation  is  an  extremely  difficult  one. 

The  feasibility  study  however  indicated  that  rms  acceleration  can  be 
correlated  with  subjective  response  of  a  ride  Jury. 
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b.  No  single  number  obtained  by  averaging  the  entire  frequency 
spectrum  provides  an  adequate  description  of  ride.  This  is  due  to  the 
number  of  different  sensing  mechanisms  involved  in  the  ride  evaluation  and 
to  the  fact  that  these  sensing  mechanisms  are  frequency  dependent;  they 
thus  "color"  the. Judgment  of  ride  quality.  A  correlation  of  the  square 

of  the  acceleration  in  each  of  the  basic  low  frequency  bands  with  subjec¬ 
tive  response  should  give  meaningful  ride  criteria  and  determine  the 
relations  between  acceleration  variance  and  ride  sensation. 

c.  It  has  also  been  concluded  that  human  response  can  be  pre¬ 
dicted  from  acceleration  variance  as  a  function  of  frequency  for  steady 
state  vibration  excited  by  stable  ground  roughness. 

Theoretical  dynamic  aspects  of  vehicle  systems 

31.  The  theoretical  dynamics  aspects  of  the  problem  are  concerned 
with  the  verification  and  development  of  mathematical  models  of  the 
vibrational  response  of  vehicles.  Although  the  modeling  methods  are 
well  established,  considerable  work  is  required  before  acceptable  agree¬ 
ment  can  be  achieved  between  actual  and  predicted  results. 

32.  The  problem  of  defining  representative  models  for  vehicle 
systems  has  been  studied  by  many  investigators  for  many  years.  Most  of 
the  publications  in  this  area  pertain  to  model  design  to  answer  specific 
problems.  The  literature  can  be  divided  into  two  different  categories: 

a.  Attempts  to  subject  models  to  artificial  inputs  and  to 
validate  the  resulting  predictioi  jy  constructing  a  similar  input  for 
the  vehicle. 

b.  Attempts  to  define  the  inputs  statistically  and  analyze  the 
behavior  of  these  models. 

33.  Most  of  the  early  analyses  pertain  to  deterministic  inputs, 
such  as  step  functions  or  sine  waves ,  and  describe  the  response  of  a 
vehicle  system  to  these  inputs.  The  method  of  analysis  used  is  either 
an  analog  computer  simulation  of  the  system  or  an  analysis  in  the  fre¬ 
quency  domain  using  transfer  functions  calculated  for  the  linear  vehicle 
system.  While  deterministic  inputs  have  yielded  direct  computer  response 
time  traces,  there  is  in  the  literature  no  available  information  on 
statistical  outputs  calculated  for  cross-country  military  vehicles.  A 
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number  of  airplane  structures  have  been  analyzed  statistically  and 

the  methods  employed  ore  of  general  interest  in  this  study.  Two  probleos 

exist  in  the  interpretation  of  previous  experiments  in  the  present  context. 

a.  Only  linear  systems  have  been  used  as  models  for  the  sta¬ 
tistical  methods  of  analysis.  Some  nonlinearities  in  the  vehicle  dynamic 
systems  are  pertinent  to  the  analysis  and  must  be  included. 

b.  All  vehicle  models  reported  in  the  literature  include  only 
a  vertical  input.  It  is  noted  from  vibration  recordings  in  an  actual 
vehicle  that  the  fore  and  aft  vibration  can  be  equal  to,  and  in  some 
cases  larger  than,  the  vertical  vibration. 

3**.  The  previous  methods  of  computer  modeling  have  been  extended 

to  allow  direct  interpretation  in  terms  of  random  vibration.  This 

extension  is  primarily  one  of  defining  the  input  and  interpreting  the 

output.  The  models  used  are  similar  to  those  in  the  literature.  Two 

different  methods  of  modeling  were  used  in  different  portions  of  the 

.study.  The  first  is  a  time  domain  analysis  which  will  allow  a  direct 

prediction  of  the  output  power  spectral  density,  and,  if  necessary  cross 

$ 

spectral  densities,  for  a  nonlinear  system.  The  second  is  a  digital 
computer  program  which  produces  the  same  results  for  a  more  detailed,  but 
linear,  system  with  the  computation  performed  in  the  frequency  domain. 

35.  A  feasibility  study  of  a  linear  vehicle  model  of  the  M3T 
vehicle  waB  chosen  for  analysis  using  both  of  the  techniques  described. 

Fig.  23  shows  a  measured  power  spectral  density  for  the  M37  vehicle  in 
cross-country  operation.  Fig.  2h  shows  the  output  acceleration  power 
spectral  density  predicted  using  the  frequency  domain  approach,  with  a 
constant-level  velocity  input.  Fig.  25  shows  the  output  power  spectral 
density,  calculated  from  the  time  domain  approach,  for  the  same  system 
with  only  one  rear  wheel  excited.  These  three  curves  indicate  that  the 
two  predictions  of  the  computer  model  output  have  essentially  three  peaks. 
The  lowest  peak  is  due  to  the  bounce  motion  of  the  vehicle,  and  the  second 
peak  is  due  to  the  combination  of  the  pitch  and  roll  motions  of  the  vehicle. 
The  third  peak  at  approximately  7.5  cps  is  a  combination  of  all  wheel  hop 
motions  of  the  vehicle,  i.e.  the  bounce  and  roll  motions  of  both  axle 
assemblies.  These  same  peaks  can  be  seen  in  fig.  23  together  with  some 
secondary  peaks  which  might  be  attributed  to  second  harmonics  resulting 


from  nonlinearities  in  the  system.  Comparison  of  these  graphB  shows 
that  the  frequency  of  the  peaks  and  the  ratio  of  amplitude  compare- 
■quite  well  for  the  two  different  computer  model  analyses  and  for  the 
measured  vehicle  response. 

36.  The  conclusions  made  on  the  basis  of  the  literature  review, 
the  development  of  the  modeling  techniques,  and  the  example  of  the  use 
of  these  techniques  are  as  follow: 

a.  Although  the  model  may  have  been  oversimplified  for  the 
analysis  used  in  the  example,  it  definitely  shows  the  feasibility  of 
using  these  modeling  techniques  for  analysis  of  vehicle  vibration  and 
prediction  of  power  spectral  density  for  interpretation  using  subjective 
criteria. 

b.  The  time  domain  approach  is  valid  for  comparison  of  general 
vehicle  configurations .  Thus,  by  treating  the  vehicle  body  as  a  rigid 
mass  and  analyzing  only  the  general  suspension  system  configuration,  we 
can  study  the  differences  between  major  vehicle  configurations.  The 
frequency  range  of  interest  is  approximately  that  from  0.5  to  10  cps 
(i.e.  to  a  frequency  above  wheel  hop  resonance).  While  this  type  of 
analysis  is  most  general  from  the  standpoint  of  including  noniinearities , 
etc.,  it  is  limited  in  degrees  of  freedom  from  the  standpoint  of  avail¬ 
able  analog  computing  facilities  in  most  installations,  and  digital 
solution  of  the  nonlinear  equations  in  the  time  domain  is  difficult  and 
time  consuming  when  the  system  becomes  very  large, 

c.  The  digital  computer  approach  for  analysis  in  the  frequency 

% 

domain  can  be  used  specifically  for  areas  of  vehicle  design  where  more 
detail  is  needed  in  the  model.  While  the  .  Jdel  must  be  linear,  it  is 
useful  for  determining  the  significance  of  attachment  points  of  various 
suspension  elements,  suspension  geometry  considerations,  etc.  The  fre¬ 
quency  range  of  interest  is  from  approximately  0.5  to  25  cps  and  this  model 
will,  in  general,  include  the  power  plant  assembly,  the  seat  on  its  mount¬ 
ing  system,  the  fundamental  bending  and  torsional  modes  of  the  vehicle 
frame  structure,  and  detail  of  the  suspension  systems.  This  type  of 
model  can  be  easily  incorporated  on  many  existing  digital  computers. 
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d.  No  solution  has  been  offered  to  the  fore  and  aft  input  to 
the  wheel,  and  it  is  3trongly  urged  that  development  be  initiated  to 
define  a  pertinent  input  of  this  nature.  It  is  not  difficult  to  include 
the  pertinent  degrees  of  freedom  in  the  computer  models  once  the  input 
is  defined. 

t 

Terrain  measurements 

37.  The  theoretical  dynamics  aspect  of  this  study  defines  the 
measurement  of  terrain  profile  necessary  for  analysis  of  vehicle  systems. 
For  the  frequency  domain  approach,  it  was  found  that  the  power  spectral 
density  of  the  profile  between  each  pair  of  wheels  was  necessary  as  an 
input.  In  the  time  domain,  a  continuous  time  trace  of  the  profile  is 
necessary;  and,  in  reality,  it  should  be  two  time  traces  separated  (in 
measurement)  by  the  tread  width  of  the  vehicle.  To  accomplish  these 
measurements,  a  brief  review  of  the  characteristics  of  the  measuring 
devices  available  was  made.  It  appears  from  the  above  criteria  that 
each  profile  should  be  measured  as  a  continuous  analog  function  and 
recorded  separately  on  two  tracks  of  a  tape  recorder,  which  would  allow 
simultaneous  playback  of  both  signals.  The  power  spectral  density  and 
cross  spectral  density  statistics  can  be  used  as  inputs  for  the  frequency 
domain  portion  of  this  analysis  and  will  permit  a  convenient  methow  of 
comparing  different  terrain  profiles  in  a  statistical  sense.  For  the 
time  domain  program,  including  nonlinearities,  it  is  necessary  to  vise 
the  actual  terrain  profile. 

38.  A  survey  of  the  existing  devices  to  measure  terrain  geometry 
was  made.  The  methods  identified  in  the  literature  can  be  divided  into 
three  general  categories:  (a)  profilometers ;  (b)  roughmeters;  and 

(c)  rod  and  transit.  The  first  two  are  automatic  mechanical-electrical 
devices,  whereas  the  last  device  is  manually  operated. 

a.  Profilometers  provide  continuous  graphic  profiles  using 
a  moving  average  reference  plane  formed  by  the  vheels  of  the  vehicle 
carrying  the  profilomtster.  The  me  asm  :rcent  is  between  this  plane  and 
a  single  profile  measuring  wheel.  The  geometry  of  the  road  wheel  pre¬ 
cludes  measurement  of  frequencies  below  a  certain  vavelength. 
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b.  Roughmeters  are  similar  to  the  profilorneter  except  that  the 
continuous  graphic  profile  is  obtained  from  a  wheel  assembly  which  acts 
as  a  seismic  device.  Ry  measuring  the  difference  in  displacement  between 
the  bed  of  the  carrying  device  and  the  wheel  assembly  forming  the  seismic 
device,  it  is  possible  to  record  continuously  the  terrain  profile. 

c.  Rod  and  transit  method  employs  conventional  surveying  tech¬ 
niques.  While  slow,  tedious,  and  expensive,  it  is  the  only  method  for 
achieving  the  necessary  level  of  accuracy  for  defining  terrain  profiles. 

39.  It  has  been  concluded  on  the  basis  of  this  study  that  the  only 
presently  available  method  of  obtaining  accurate  terrain  profile  data 
suitable  for  analysis  of  vehicle  systems  is  the  slow  and  expensive  rod 

and  transit  method.  It  is  suggested  that  a  device  acceptable  for  measuring 
terrain  profiles  automatically  should  be  capable  of  measuring  ground  wave¬ 
lengths  between  0.011  and  2.73  cycles  per  foot  as  a  continuous  analog 
recording  of  two  tracks. 

Conclusions 

40.  As  might  have  been  expected,  the  two  approaches  to  the  problem 
of  describing  and  defining  the  dynamic  response  characteristics  have  not 
led  to  a  clear  judgment  of  the  superiority  of  one  over  the  other.  They 
have,  however,  forced  the  schools  of  thought  to  shift  their  ground  slightly. 
The  FMC  approach,  which  traces  force  and  motion  in  real  time  through  the 
vehicle  suspension  system,  has  been  demonstrated  to  be  a  useful  design 
tool.  In  fact,  the  FMC  Corporation  is  currently  using  it  to  design  the 
suspension  system  of  its  new  armored  personnel  carrier,  and  it  appears 
safe  to  predict  that  the  new  machine  will  be  a  substantial  improvement 
over  the  M113.  The  force-motion  system  does  not,  of  course,  yield  a 
description  of  ride  which  can  be  readily  correlated  with  human  response. 

For  this,  the  Chrysler  Corporation  statistical  system  appears  to  be  the 
only  solution  in  sight. 

41.  Thus,  it  appears  that  the  dual  approach  has  revealed,  among 
other  things,  that  there  are  two  separate  and  distinct  problems  instead 
of  one  as  had  been  previously  assumed.  One  problem  is  to  define  force 
and  motion  in  a  suspension  system  in  such  a  way  that  the  relations  can 
be  used  to  improve  the  design  of  suspension  systems.  This  the  FMC 


Corporation  mathematical  model  appears  to  do.  The  other  problem  is  to 
describe  motion  in  objective  terms  in  such  a  wny  that  the  mathematically 
derived  terras  correlate  directly  with  human  response.  While  a  satis¬ 
factory  solution  to  this  problem  has  not  yet  been  achieved,  it  appears 
that  the  Chrysler  Corporation  statistical  procedure  will  be  suitable 
for  this  purpose. 
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Table  1 

Deflections  and  Deflection  Ratea  Used  in  Vehicle  Model 


TIRE  DEFLECTIONS: 


SUSPENSION  DEFLECTIONS: 


s,-xx,- 
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a*-  x+Jifd-*-^  -7*  A* 
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Table  2 

Equations  of  Motion  for  4-Wheel  Vehicle 

BODY  BOUNCE  EQUATION: 

(-K,-Kt)X„+(-C,-CI)X,i*(-Ks-K0XM  +  (-C,-t4UM 

+(C»-C4)^,4+(-X|-Xi-K,-K4)X+(-C1-Ct-CrC4)t-j? 
♦[(-K.-KOJ-ff  ♦(K4*K4)^^r]e-*[K-c,)ll(4(cltC4UJr]e 
*^,-Ki^*C-X»+K4)^]y+[(C,-<l^*(-Ci+C4^]v  *  0 

BODY  PITCH  EQUATION: 

C-K,-K i)M(  Xtt  +(-C,-Ci)Uf  X„  +  fKj+K4 )/Jr  X,4 

♦(•K!«4)^4^44(-C,*<4)«r^M 

*  ^  K,- Kt)l  X(  +(K,*K4M  tr  j  X  ♦  [(-  C ,  -  Ct)  X  i4  +(£,  ♦£«)«,]  % 

♦[(-K1-XOUirtl*(-Ks-K4Xiir)']e*[(-C,-C,')(/l>,)S(.CrC,)«/r)']i 
-ly 8  *[(K, -X0  £  ♦  (Ki-K^JtJ,  i.]  y 

♦  [(C,-CiVit^4(C%.C4^JUr^]lf  =0 
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Table  2  (con.) 


BODY  ROLL  EQUATION: 

+(K|+K0(^r)  0!t+  (C,+C0(y) 

♦[(t.-C,)^  +  (-C4*C4^]{  +  [(K1-Kl^Wt*(VK4^iA]e 

+[(<V<0^.*V(c*-M  9+[(-X,-K1#)+(-VK4#)I1y 

■*[(-c -Ct^V-Ci-^yj'f-XiY  -  0 


FRONT  DIFFERENTIAL  BOUNCE  EQUATION: 

KKlXXl*CC,xici^KKtXXt^CC1^t«*-(-KKl-KKt-Kl-Kl)X,t 

•K-CCrCCt-C1-Cj)X,t-^1i4[(-KK,tKK,)|+(-K,*Ki^]^t 

+[c-cc,*cct>| Jf(-c,-Ci)£ 

+  (-K,-K t)JUf  8*(-C,-Ct)  Jfif  6  +(K,-Ki> ^  V 

FRONT  DIFFERENTIAL  ROLL  EQUATION: 

KK,  j  XX,  +  «,  |  xV, -KKti-XXt -CCt ^-XX, 
4-[f-KK,+KH,>|--<-K)*K,V^rj<.t  +  |'a1*CC,)|*(-C1*Cl')^]xlt 

♦[{-KXrKKl^)V(S(,-K1](ii)2j^,l«[(.«rCC,||)V(-Crf^),]?-t 

~(Ia)l  ^il  +(-K,-*.Kt)  ^  X  ■*'("£, +  0^)^*  X*(-Kt4Xv)Y  ^4  ® 

®  *(x,*K,)(^y  y*(c,*cv 
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Table  2  (con.) 


REAR  DIFFERENTIAL  BOUNCE  EQUATION: 

KK4XX44CC|XX5+KK4XX4  +  (.C4X\4  «‘(-KK*-KK4-Ki~K4)X|f 
4  ;-<Cj-CC4-CrC4)XH- ^i<j4+[(KKi-KK4)|«{lCrK4.')^]^4 

+[(cc»-co|*(c,-«4>  £  j  +(-K,-K4)j!+(-c»-<^f  4-(K»*K  OM,e 

’(Cj»C4)^,l+(-K,*K4)Y:V'  +  (-C4*C4')i.'f  .0 

REAR  DIFFERENTIAL  ROLL  EQUATION: 

-KK.,|xX,-CCj|  X*Xj  +KK4|xX4  +  CC4ixV4 
♦[(KKi-KK^I  +(XS-K4)  *f}  XM*  [(CC,-CC4)f  +(C,-C4)  lf]  X,4 
♦  [(-XXj-KK4^)V(-K»-K4^)*-]^4  -  (lu)t'l?j4 

*[(-«1-CC4^)‘*(-C,-C4l(t)t]?l«A  +(«*-**) 

■*•  (C»*C4^  i  +(,-X»+K4)  f  Hr  6  +(-  C»*C4>  £  Ur  8 

+  (KS+K4)(^-)*V'  +  (C44C4)(4:)^  -0 
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Body  Bounce 


Body  Pitch 
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Fig.  1.  Method  of  extending  surface  geometry  profiles  for  Fourier  analysis 


Fig,  3.  Site  152*  rice  field  bund 


Fig.  1*.  Site  152 •  Field  prepared  for  row  crop 
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THREE  DIMENSIONAL  MATHEMATICAL  MODEL  ©F  MVf  TYPE  TRUCK 
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Fig.  15.  Schematic  of  general  model  of  wheeled  vehicle 
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Fig.  19.  Sample  of  course  constructed  for  ride  dynamics  model 
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Mean  square  acceleration  per  cps  (g2/cps  x  103) 


Mean  square  acceleration/ cps 


Mean  square  acceleration/cps  (in./sfc:;2)2  per  CpS 


*  l  I  I  i  I  J  -  l—  i.  .  J _ I _ L-  1—3  . 

(Predicted  at  a  rigid  vehicle  seat 
position  using  frequency  domain 
—  approach)  \&C? 
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Fig.  2k.  Output  acceleration  power  spectral  density  for 

constant  level  velocity  input  power  spectral  density 

P-  (acc)  =  u)2  P_  (vel) 


square  acceleration/cps  (g2/cps  x  103) 


Appendix  A 


ELECTRONIC  COMPUTER  PROGRAM  ABSTRACT 


TITLE  OF  PROGRAM 

Fourier  Series  Curve  Fit  for  Unequally  Spaced  Points 


PREPARINO  AGENCY 


AUTHOR  IS) 

status  of  program 

J.  F.  Smith 

March  1965 

RH  ASK 

MOD 

STATIC 

OP 

To  compute  the  amplitude  coefficients,  c1(  and  the  phase  angles,  o^,  for  N 
terms  of  a  Fourier  Series  to  fit  a  set  of  M  equally  or  unequally  spaced 
points  with  linear  interpolation  between  points.  The  output  is  used  in  an 
equation  of  the  form  M 


+  I  ci  sin  +  ai ) 


B.  PROGRAM  SPECIFICATIONS 

The  program  will  accept  up  to  a  maximum  of  100  coordinate  pairs  as  input  and 
it  will  compute  up  to  52  terms  in  the  Fourier  Series.  A  control  card  contains 
the  number  of  coordinate  pairs,  the  number  of  series  terms,  and  flags  to 
obtain  output  options . 

The  program  is  written  in  FORTRAN- I I  and  uses  184110  words  of  storage 
excluding  subroutines. 


C.  METHODS 

The  Euler  formulas  are  computed  using  linear  interpolation  between  points  to 
obtain  the  aj^  and  b*  coefficients.  Then  we  have 

ci  c  <ai2  +  bi2> 


where  1  s  1 ,  2 ,  . . . ,  N. 

(See  reference  for  derivation  of  this  equation) 


D.  EQUIPMENT  DETAILS 

GE  225  computer,  8K  memory,  card  reader,  punch,  and  on-line  printer. 


E.  INPUT*  OUTPUT 

Coordinate  pairs  are  punched  in  9-column  fields,  4  pairs  per  card;  l.e.,  x^, 

y  ,  x  ,  y  ,  x  ,  y  ,  x.,  y. ,  etc.  The  data  with  a  decimal  point  can  be 
X  2  2  3  3  4  4 

anywhere  in  the  field.  (See  write-up  for  complete  input  format) 

Output  is  printed  listing  with  optional  punch  cards  of  either  the  normal  Euler 
coefficients  and/or  the  amplitude-phase  angle  coefficients. 


F.  ADDITIONAL  REMARKS 

Original  code  was  written  by  A.  F.  Zettlemoyer  and  C.  W.  Bowman,  Manufacturing 
Research  Space  Guidance  Center,  Oswego,  New  York. 

Reference  fur  the  equation  derivation  is  C.  R.  Wylie,  Jr.,  Advanced  Engineer¬ 
ing  Mathematics ,  McGraw-Hill,  1951,  pp.  113-130. 

Write-up,  listing,  source  deck,  and  binary  deck  available. 


ENG  FORM  aaa« 
t  AUO  CO 


m 

k 

lv- 

m 

IV* 

ki¬ 

rn 


few  mo'Wifo) 


AAKVIOUS  KOITIONS  ARK  OSIOLXTC. 
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FOURIER  SERIES  CURVE  FIT  FOR  UNEQUALLY  SPACED  POINTS 


Specifications 

Coded  by 

Entry 

Type 

Storage 

Purpose 


J.  F.  Smith,  March  1,  1965 
Main  Program 
FORTRAN- II 
184110  (GE-225) 


Given  .a  series  o£  coordinates  representing  a  function,  to  compute  the 
Euler  coefficients  for  a  given  number  of  terms  of  a  Fourier  Series  that 
best  represent  the  function. 


Use 

1.  The  function  to  be  fitted  must  be  expressed  in  a  series  of  coordinate 
pairs  (up  to  100).  These  pairs  may  or  may  not  be  equally  spaced. 

2.  An  input  value,  N,  indicating  the  number  of  terms  to  be  computed  (up 
to  52)  aud  a  set  of  flags  controlling  output-  options  are  entered. 

3.  The  output  is  a  listing  of  the  coefficients  and  the  calculated  "fit" 
of  the  input  data.  Optional  punch  card  output  is  also  available. 


Method 


Straight  lines  are  calculated  between  each  set  of  data  points  in  the 
slope,  c,  and  intercept,  cc,  form.  Each  of  these  linear  equations  is  used 
as  the  f(x)  in  the  interval  for  the  calculation  of  the  Euler  coefficients. 

The  following  integrations  are  performed  for  a  ,  the  constant  term;  a^, 
the  cosine  coefficient;  and  b^,  the  sine  coefflclen?. 


a 


o 


f  (.c)  dx 


a 


i 


f(x)  cos  (~  x)  dx 


(1) 


(2) 


(3) 


/P 

f(x)  sin  (~  x)  dx 

U 

where  2p  is  the  fundamental  period  (2p  = T) ,  f(x)  =  cx  +  cc,  and 
i  a  1,  2,  ....  N. 


To  speed  up  evaluation  when  using  the  coefficients  to  evaluate  a  forcing 
function  in  a  dynamic  simulation,  an  alternate  form  of  Fourier  Series  is  used. 

a  N 

*(t)  =  y  +  Z  ci sin  t‘+  aiJ  (4) 

1=1  T 


where  t  =  independent  variable 
f(t)  =  dependent  variable 

a0  =  twice  the  average  value  of  f(t) 
c*  =  amplitude  coefficient  for  the  ith  term 
O'!  =  phase  angle  for  the  i^  term 
T  =  the  fundamental  period 
N  =  number  of  terms  in  the  series 


To  compute  and 


for  use  in  this  equation,  we  have 


c 


l 


2.i 

+  b^  ) 


(5) 


(6) 


where  a^  and  b^  are  from  equations  (2)  and  (3)  and  i  =  1,  2,  . K. 


Data  preparation. 


Card  1 


?1  Hollerith  characters,  cols  2-72 


Card  2 

Cols  3-10 

M 

11-20 

FC* 

?l-30 

CARD 

33-40 

ABLIST 

41-50 

ABOARD 

Card  3,  3+n 

Cols  3-9 

T(D 

10-18 

AI(l) 

19-27 

T(2) 

28-36 

A)  (2) 

37-45 

T{3) 

46-54 

AI  (3) 

55-63  • 

T('l) 

64-72 

Af  (4) 

A-3 

(rt.  adjusted  in  the  field) 
(anywhere  in  the  field  with 
a  decimal  point) 


(anywhere  in  the  field  with 
a  decimal  point) 


Continuo  T(i)  ,  AI(i)  values, 
4  pairs  per  card  for  K  pairs 


FORTRAN  Symbol  Pefinlt ion 

M  manor*.-  of,  pairs  of  coordinates 

?•"&  number  of  terms  in  computed  series 

CARO  1.0  indicates  punch  ampl.  coeff,  and  angles 

0.0  indicates  no  punch 

ABLIST  1.0  indicates  print  a^  and  bj 
0.0  Indicates  no  print 

ABOARD  1.0  indicates  punch  a^  and  b^ 

0.0  indicates  no  punch 

T(i)  time  (or  distance)  coordinate,  i  =  I,  M 

AI (1)  amplitude  (or  elevation)  coordinate,  i  =  1,  M 


Output 

The  71  Hollerith  characters  are  printed  in  the  heading.  The 
amplitude  coefficients  and  phase  angles  are  printed  and  an  option  allows 
a^  and  b^  to  be  printed. 

The  tabulation  shows  the  input  data,  x  and  f(x),  with  the 
calculated  f(x)  and  error  with  an  absolute  percent.  The  full  scale 
error  is  the  maximum  absolute  error  divided  by  the  maximum  f(x). 

Punch  card  output  is  optional.  Set  CARO  =1.0  and  the  amplitude 
coefficients  and  phase  angles  are  punched  out  in  groups,  4  per  card,  with 
E16.7.  The  output  looks  as  follows; 


Card 

1 

Cols 

1-16 

ao 

bb±0 .FFFFFFFE ±xx 

17-32 

c(l) 

33-48 

c(2) 

. 

49-64 

c(3) 

Card 

2+n 

Cols 

1-16 

c(4) 

Continue  c(i) ,  4 

per  card, 

17-32 

c(5) 

for  N  values . 

33-48 

c(6> 

49-64 

c(7) 

Card 

3 

Col  8 

1-16 

T 

bb±o . FFFFFFFE ±xx 

17-32 

«(1) 

33-48 

or(2) 

49-64 

»(3) 

Card 

4+n 

Cols 

1-16 

«(4) 

Continue  or(i).,  4 

per  card, 

17-32 

<*(5) 

for  N  values . 

33-48 

«(6) 

49-64 

“(7) 

A-h 


Tho  optional  punch  output  of  and  is  the  same  as  above ;  aQ  and 

T  are  punched  and  b  =  0.  Set  ABOARD  =1.0  and  output  Is  as  follows: 

n 


Card  1 

Cols  1-16 

ao 

bb±0 . FFFFFFFE  ±xx 

17-32 

a(l) 

33-48 

a  (2) 

' 

49-64 

a  (3) 

Card  2+n 

Cols  1-16 

a  (4) 

Continue  a(i) ,  4  per  card. 

17-32 

a  (5) 

for  N  values . ’ 

33-48 

a  (6) 

• 

49-64 

a  (7) 

Card  3 

Cols  1-16 

T 

bb±0 .FFFFFFFEiXX 

17-32 

b(l) 

33-48 

b(2) 

49-64 

b(3) 

Card  4+n 

Cols  1-16 

b(4) 

Continue  b(i) ,  4  per  card, 

17-32 

b<5) 

for  N  values. 

33-48 

b<6) 

49-64 

b<7) 

GE-225  User  Kotes 

1.  If  function  being  fitted  is  to  be  used  in  a  cyclic  manner,  then  the 
first  and  last  pairs  of  data  points  should  agree  in  amplitude  to  avoid  a 
discontinuity. 

2.  The  statements  to  read  the  coefficients  punched  by  this  program  into 
another  FORTRAN- 1 1  program  knowing  N  are : 

10  FORMAT  (4E16.7) 

READ  10,  AO,  <C(I),  I  =  1,  N) 

READ  10,  T,  (ALPHA (I)  ,  1*1,10 

OR 

READ  10,  AO,  (A (I) ,  I  =  1,  N) 

READ  10,  T,  (B(I),  I  =  1,  N) 


References 

This  program  was  originally  written  by  A.  F.  Zettlemoyer  and  C.  W.  Bowman, 
Manufacturing  Research  Space  Guidance  Center,  Oswego,  New  York. 

The  derivation  of  the  alternate  form  of  the  Fourier  Series  may  be  found  in 
C.  R.  Wylie,  Jr.,  Advanced  Engineering  Mathematics,  McGraw-Hill,  1951,  pp. 
113-130. 
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FOUR  I ER  SERIES  CURVE  FITTING 
UP  TO  1GU  UNEUUALLY  SPACED  COORDINATE  PAIRS  ALLOWED  AS  INPUT 
UP  TO  52  TERMS  COMPUTED  AS  OUTPUT 

output  print  of  amplitude  coefficients  and  phase  angles 

OPTIONAL  CARD  OUTPUT'oF  AMPLITUDE  COEFFICIENTS  AND  PHASE  ANGLES 

OPTIONAL  PRINT  OF  THE  SINE  AND  COSINE  COEFFICIENTS 

OPTIONAL  CARO  OUTPUT  OF  THE  SINE  AND  COSINE  COEFFICIENTS 

BY  A.K.  ZETTLEMOYER  AND  C.W.  BOWMAN#  MANUFACTURING  RESEARCH 
SPACE  GUIDANCE  CENTER#  OSWEGO#  NEW  YORK 

DIMENSION  Tl 100 1 #AI 11001 ,AI52] ,AA 152) #CFl53) *  All  531 
EO-JIVA!  ENCF  ICFIU.U),  IAL11J#PTT1 

105  FORMAT [ 110, bE10.8) _ _  _ _ _ _ 

1 jo  KURMAT  I6F9.4] 

3o7  FORMAT  ( **E16  •  7  ]  .  -  • 

l2u  FORMAT  [12X1HX16XRHFI  XU1X8HCAL  F  I  X  )  12X5HERROR10X7HPERCENT  ] 

121  FORMAT  (///  41X  10HT  ABULAT I  ON  /}.,  . .  . . 

3.22  FORMAT  flh  ,  1P4E18 . 6, 2PF11 .  ?  ] 

3  23  F0RMAT(1H1.29X#.28HF0.uRIER  SERIES.  CJ.IR.V.E..F.1.TTI.NG..././/J . . 

127  FORMAT I72rt 

1  J.  - . .  -  . . 

32o  FORMAT [ 1H0 ] 

1002  FORMAT  (E20.7#E18.  / » 41X  #  15)  . . .  . 

1  j  u5  FORMAT,"  1?X4HA[0]15X1HT1 

M-u7  FORM  A  T  I  /8X12HAMPL  •_  _C_06FF  .  6X1 1H  ANGLE  ...I  RAD  1.4  5X1 HN  1.  ..  . . . . .  . . 

ICON  FORMAT  (  iliO  #  6X ,  19HFULL  SCALE  ERROR  =  ,2PF5.21 

3JG9  FORMAT I/8X12HAMPL.  COEFF . 6X11HANGLE  f RAD ) 10X4HA l N 1 15X4HB l N 1 12X1HN ) 
1010  FORMAT IE20.7#3E16.7*5X# 151 

INPUT  COMPUTATION  LIMITS  AND  FLAGS'" 

■»  » 

30  READ  127 

RE  AO  105#M#FCN#CARO#  ABLIST,  ABOARD . . . . . 

N  a  M  -  1 

NMAX  =  FCN  ....  .  .  . . .  . . 

;  IN^UT  COORDINATE..  UA.T.A _ ... _ 


READ  106#  l  T  (  I  1  #AI  II.)  #  I=1»M] 

PTT=TIM]-Tfll 

PT=.5*PTT 

DO  701  1=1.52 

A  l  I  ]=U. _ _ _  . 

7ul  A  A ( I ) =  0 . 

XP!=3. 14159265 

FN  =  1. 

U=0.  . 


COMPUTE . SINE _ANp_COSLNE_COEFF ICI ENTS _ 


00  45  ,J  =  1»NMAX 
L  =  2 

<#j  43  I = 1 #  M 


no  o  o  o  o  o  o  non 


CslAI ILJ-AI I  1 1 1 / l T I L ) -T  t  I  10 

CC  =  AI  I  I  J  -C*T  (  I  I  . .  .  .  .  . . . 

IFIFN-1..0)  22,  22,  23 

22  Ssl.O/PT  *1C*TILI*TILJ*0.5+CC*T|L|-c*TC 11*TII!*0.5). 
Vsl.O/PT  * I-CC*T K 1 1 J 

UsU  +  S  +  V  _ _ _ _ _ 

23  XF=FN*XPI 

XHi)=OPT/  (XF*XF  ) _ _ , _ 

XRjsC/XF 

CCX=CC/XF  _ 

X  F I s  X  F  *  T { IJ/PT 

XFLsXF*T(L]/PT  _  _ 

SXSsSIwF(XFL) 

SX3  =  SINF(XFIJ _ _  _ _  _ ...... 

osnssxe-sxg* 

CXEsCORFIXFLl 
CXGsCOsF ( XF I  1 
OCN=CXE-CXG 

A l JJsXHD*DCN+XHI*iTlL J*SXE«T t 13 *SX6 ) +CCX*DSN+A t J 1 
AAt  J]=XHD*  D  S  N  -  X  HIM  T  JJL.l.  *  C  X  E  -  T  I. 1  j.*  J -_C  C X  *  fl  C  N + A  A  .J.  JJ . 

L*L  +  1 

41  CONTINUE  .  _ _ _  _ _ _ _ _ ... 

FNsFN+l . 

4t>  CONTINUE  _ _ _ _ _ 

PRINT  HEADINGS _ 

PRINT  123 _ _ _ 

PRINT  127 

PRINT  12*  _  _  _ _ 

PRINT  1005 

...PRINT.  10.02,  U*  PIT. _ _ 

COMPUTE  AMPLITUDE.  COEFFICIENTS.  AND..PHASE  ANGLES 

DO  1040  K=1,NMAX  _  _ 

COAMP  =  SORTF  l  A  IK)  ★  A  t  K  J  ♦  AA  (K)  *A~A  IK)) 

AlIK+11=ATANF|A(K1 /AA l K ) ) _  _ 

1037  IF  I  A A { K ] 11038,103971039 
1IM6  CO AMP  =  -COAMP 
1.13 9  CF  ( <  +  l )  sCOAMP 

11)40  CONTINUE  . .  .  _ _ _ _ _ 

PRINT  OUTPUT _ J. _ 

AA  t  NMAX 1 =  0  ,  G 

IFIABLISTU041, 1041, 1042 

1041  PRINT  1007 

PRINT  1002, tCFIK+l) , ALIK+1] # K # K«1 # NMAX ] 

'  GO  TO  1043  _  _  _ _ 

1042  PRINT  1009 

PRINT  1  01  U,  (CF  (K  +  l  ]  ,  ALJK  +  1]  ,  A  l  K  AA  IK  ]  ,  K  ,  K»1 ,  NMAX  ) 

1043  CONTINUE 

PUNCH  CARD  OUTPUT  . . 

p 

IF  I  CARD  ]  l’G45»'l(M5,  in46 
1f;4ft  NCAH0  =  NmAX>1 

PUNCH  107, ICFIK] #K=1,NCARD) 

PUNCH  107#  (AUK)  ,Ksl,NCARD]  ..  .  . . . . 

A- 7 


o  o  o 


:J45  CONTINUE 

lMARCA«n)1047, 1047, 104b  _ _ 

104t)  PUNCH  107,u/tA(K’i/KVt,NMAXV 

PUNCH  10?,PTT,  tAAlK),K=l,NMAXl 
1047  CONTINUE 

COMPUTE  AND  PRINT"  FIT ‘anO  ERROR  OF  F if 
XI  =  ,5*lj 

XDIFFsn.0..„  _ _ _ _ 

AMAX=U . 0 

1)0  1060  1*1, M  _  _•  _ 

X  =  X1 

DO  1050  K*1,NMAX _ _ 

Z  =  K 

CO  AMP  =  CF  [  K*1 1  ....  .  _ _ 

YbALIK+1) 

1051)  XsXfCOAMP*SINF  IZ*XPI*T  (  I  1/PT+Yl  _  _ 

1 F II  -  21  805,  806,  806 

805  PRINT  121 _ _ 

PRINT  120 

806  DIFF  =  X  “  AIII1  _ 

IF ( A I l  1  1  1808,807,808 

e07  PCD  I F  =  0 . 0  _  _ _ 

GO  TO  DO 9 

608  PCDIF  =  ABSr  IDIFF/AIJ  1 1)  _  _ 

809  PRINT  122,  T  t  I  1 ,  A  11  I  J  ,~X ,  D  IFF  ,  PCD  I F 

IF l XD IFF-ARSF IDIr FJ ] 1030,1031# 1031  _ 

1030  XDlFFsABSPlDIFF) 

1031  IFIAMAX-AHSFIAI  t  IJ_1  11032, 1035,_1035 _ 

1.(32  AMAX*ABSF  IAI  tlli 

1  :>  3  5__C  ONTINUE _ 

1060  CONTINUE 

XDIFF  =  XDIFF/AMAX  __  _  __  _ 

PRINT  1 0  Ob , XD I FF 


(M'CM 


.._  .0  JUNE  1965  J~  SMITH _ DATA  1 

1.  DATA  2 


-2.0 

.03 

_ -2.0 

..  DATA  3 

-2.0 

.07 

-2.0 

DATA  4 

_  -2. 0 

.11 

-2.0 

..DATA  5 

-2.0 

.15 

.-2.0 

DATA  6 

-2.0 

.,..19. 

-2.0 

-DATA  7 

-1.9 

.23 

•  -1.35 

DATA  8 

....  1.3 

.27 

.  ..  1.9  . 

.DATA  9 

3.11 

.31 

3.15 

DATA  10 

.  2.8  _ 

..35 

2.78 

.DATA  11 

2. 78 

.39 

2.78 

DATA  12 

2.73 

..^3 

.  .  2.70 

.-D..A.TA...13 

2.50 

.47 

2.40 

DATA  14 

_ _ .-2.0 _ 

...51 

,.-2.08. 

DATA  15 

-2.0 

.55 

-1.92 

DATA  16 

-1.95 

.59. 

-2.0 

.  DATA.  17. 

-2.0 

.63 

-2.0 

DATA  16 

-2.0 

.♦A7_ 

-2.0 _ D.A.IA  19.- 

-2.0 

DATA  20 

FOURIER' SERIES  CURVE  FITTING" . 


‘'ll 


I.RIEK  FIT  TEST  DATa':  SAMPLE  PROBLEM  A  8  JUNE' i965'  '  J?  SMITH 


)  4(0) 

-IK  76428556*00" 

|  .AMPL .  COEFF. 
j  -0.263O392E+01 
U.1266593F+01 
r  0.1619213E-01 
-0.58033316+00" 
0 .42237136+00 
-0 .74261346-01 
-0.2387C62F+Q0 
0.1519472F+00 
0.1174537F+00 
-0.1194020F+00' 
0.3952108E-01 
U. 93607756-01 
-0.7324159F-01 
-G.1998C52E-01 
f.6323731E-01_ 
-0. 35898526-01 
-0.21264906-01 
U. 37736596-01 
-U.12U0203E-01 
-u.164A09!5E-01 
]  0.16U557F-01 
O.35945'O'8E-02“ 
j-0. 126(15956-01 
|  0.1167912E-02" 
j -G . 12387996-01 


"O'.  700000  O'6  +  OO" 

“  ANGLE  (RAD!  “ 
0 • 14335466+01 
0 .13012546+01" 
0.58169756+00 
0797874326+00" 
0 . 72789406+00 
•0.10779316+01 
0.60396926+00 
•0 . 12629376+00 
0.1244368E+01_ 
0.4  5  841 3  8  6+0  0 
■0.1063192E+01 
0 . 6097020E+00" 
0.20669286+00 
0'.'i0973llE+01 
0.15510916+00 
•0*24  i’O  7  256+00 
0.1139477E+00 
•0  .’32308706+0  0' 
•0.6723846E+00 
-0.5602277E+00' 
•0^,83219666+00 
"0 .772 1 9  4  5  6 +0  0~" 
•Oj.  11732556+01 
•079*6045726+0  0" 
0.14272336+01 


A  INI' 

-0.2614571E+01 
0 .12208596+01 
0.8896656E-02 
-  0.48155856  +  0  0~ 
0.2810031E+00 
0.65422816-01 
-0.13556466+00 
-0.19139016-01 
_0. 11144096+00 
••0 . 5283850E“dl~ 
-0.34537926-01 
0.53601976-01' 
-0.15030956-01 
-0.17782356-01 
_0.9769400E-02 
'0.'857O564E-O2~ 
-0.2417847E-02 
"-0.11981196-01“ 
0.74755076-02 
0,87470066-02" 
^0.11916046-01 
-0.25079  2  i  6  -  0  2~ 
0.11622886-01 
-0.95705006-03" 
-0.1226Q55E-01 


B(N) 

_ -0.36112006+00 

0.33728156+00 
0.1352903E-01 
-0.32386416+00“ 
0.3153328E+00 
"'  -0.3513690E-01 
_^0. 19647626+00 
0.15073716+00 
0 » 3772702E-01 
-0.1070744E+00" 

_ 0.19210606-01 

0.7674138E-01 

_ -0 . 7168264E-01 

-0.91109356-02 

_ 0.6247812E-01 

-0 .34 860426-0 l" 
-0.21127006-01 
0.35784096-01" 

_ -0.93896556-02 

-0.13944636-01" 

0.10849866-01 

0.25750386-02- 

_ -0.48804306-02 

0.66936356-03 

O.OOOOOOOE^OO 


o.ononcoE-oi 

1.0000006-02 

2 . 0000006-02 
3. 00 0  0  COE-02 
j  4.00C0006-02 
i  5 , 00QO00E-02 
j  o . onoooofc-o2 
i  /.  00000 uE-02 
■  fc.  oo  11  no 06- 02" 
j  9.0000006-02 
!  1 . 0  0  0  n  0  0  fc  -  0 1 
j  l.lOOOCOE-Ol 
1.2C0100E-01 
'  1 ,3(il>00l'fc-01 
!  : .4ocoooE-of 
;  :.5uooooe-oi 
1  .  h00O00E-0l’ 
j  1.7008006-01 

l.aoonooE-oi' 

1.900000E-01 


FIX) 

-2. 000000E+ 
-2.0000006+ 

~270“00000E4 
-2.0000006+ 
-2.000000E+ 
-2.0000006+ 
-2.0000006+ 
“2.0000006- 
-270000'00fc  + 
-2.00000064 
-2.0000006+ 
-2.000000E4 
-2.0000006+ 
_^2,JOOUOOOE  + 
-  2 . 0  0  0  0  0  0  E  ♦ 
-2.000000E4 
-2.000000E+ 
-2^0000006+ 
"-2.000000E+ 
-2.0000006- 


_ TABULATION _ 

CAL  FIX) 

00  ' -2. 007796E+00 

00  -1.996401E+00. 

O’O  -i.'99'7043E  +  0'0' 

00 _  -2.0075226+00 

00  -i. 9929426+00 

00  _  -2. 001826E  +  00_ 

00'  -2.004757E+00 

00  -1.991695E+U0 

00  -2.0  063246+0  0“ 

00  _  -2.000084E+00 

00  -1.993188E+00 

00 _  -2.009312E+00 

00  -1.994437E+00' 

00  -1.997574E+00 

O'O  -2.0096326+00“ 

00 _  -1. 9889676*00 

00r'  '  -2.004853E  +  00' 

00 _  -2. 0059326+00 

00'  -1.9851156+00 

00  A-10  -2.0155296+00 


ERROR 

-7.7955576-03 

3.599301E-03 

‘2.9573946-03* 
-7.522304E-03 
7.0584176-03 
-1.825713E-03 
-4.7569576-03 
8.304773E-03 
^6,3243326-03 
-8.354709E-05 
6.812410E-03 
-9,3116466-03 
5.563444E-03 
2,4258596-03 
•9. 6322606-03' 
1.1032616-02 
-4.8533796-03' 
-5.932171E-03 
"1.4885466-02 
-1.552910E-02 


PERCENT 

"'0.39 

0.18 


4.  •  U  • '  U  '  »)  U  C  -  U 

-k.uuuuuue^uu 

-x. V951246+UU 

4.8//3266-03 

0,24 

^.luonoofc-oi 

-2.0000006+00 

-1.9793386+00 

2.0662096-02 

1.03 

k .  kOl-nOOE-Ol 

-1.9000006+00 

-1.8973636+00 

2.637347E-03 

0.14 

k .  3h(' 9  00E-01 

-1.350000E+00 

-1.224547E+00 

li254531E-01 

9.29 

k.4fti«OnE-01 

2.000000E-02 

-8.8909686-02 

-i»  089097E-01 

544.55 

k.50i1000E-01 

7.500000E-01 

7 . 649811E-01 

1.498110E-02 

2.00 

k .  6  0  0-1 0  0  6  -  0 1 

i.300000E+0d 

1. 3016866+00 

1 . 688389E-03 

0.13 

2.7000-006-01 

1. 9000006+00 

1. 8751476+00 

-2.485264E-02 

1.31 

k.dOfinoOE-Oi 

2.4300006+00 

2.467391E+00 

3.739059E-02 

1.54 

S.900n0UE-01 

3. OOOOOOE+OO 

2.9243176+00 

-7.568309E-02 

2.52 

3. 0000006-01 

3.1100006+00  ’ 

3.1321006+00 

2.210039E-02 

0.71 

3.  1000006-01 

3. 150000E+00 

3.1216206+00 

-2.837976E-02 

0.90 

•3.200000E-01 

3.0300006+00 

3.0331956+Od 

3.194928E-03 

0.11 

3.3000POE-01 

2. 900000E+00 

2.9O5005E+00 

5. 005337E-03 

0.17 

3.400000E-01 

2. 8000006+00 

2.8029906+00 

2.989691E-03 

0.11 

3.‘5n0000E-01 

2.7800006+00 

2.7843856+00 

4 . 384723E-03 

0.16 

3.600000E-01 

2. 7800006+00 

2.7794486+00 

-5.520657E-04 

0.02 

3.700000E-01 

2.780000E+C0 

2.7778866+00 

-2.114005E-03 

0.08 

3.800POQE-01 

2.7800006+00 

2.7841966+00 

4.195947E-03 

0.15 

3.9000006-01 

2.780000E+00 

2.7741296+00 

-5 « 870745E-03 

0-.  21 

■4  .  ooonooE-oi 

2.7700006+00 

2.773851E+00 

3.850516E-03 

0.14 

a.ioooooe-oi 

2.760000E+00 

2.758091E+00 

-1.909409E-03 

0.07 

4.20onoofc-Gi 

2. 730000E+00 

2.7251496+00 

-4.850894E-03 

0,18 

4.3000006-01 

2.7000006+00 

2.7075036+00 

7.50285CE-03 

0.26 

4.4(iUO0nE-01 

2.6600006+00 

2.6529566+00 

-7.044449E-03 

0.26 

4 .50(1-1006-01 

2.600000E+00 

2.5882636+00 

-1.173666E-02 

0.45 

4.6000006-01 

2.5000006+00 

2.537562E+00 

3,7561676-02 

1.50 

4.7000CC6-01 

2.4300006+00 

2.2381466+00 

-1.618540E-01 

6.74 

4.8000006-01 

1.2000006+00 

1.2747986+00 

7.479822E-02 

6.23 

4.90000CE-01 

-2.0000006-01 

-3.446057E-01 

-1 . 446057E-01 

72.30 

5.000000E-01 

-2.000000E+00 

-1.7434406+00 

2.565604E-01 

"12.83 

5.1000006-01 

-2.080000E+00 

-2.1526856+00 

-7.268480E-02 

3.49 

5.2000006-01 

-2.0900006+00 

-2.049929E+00 

4 . 0  07091E-02 

1.92 

5.3000006-01 

-2.050000E+00 

-2.054098E+00 

-4.097831E-03 

0.20 

5.4000006-01 

-2.0000006+00 

-2.010779E+00 

-1.077857E-02 

0.54 

i  5 . 5 1 0  0  0  0  E  -  0 1 

-1.9200006+00 

-1.907840E+00 

1.216040E-02 

0,63 

JS.frOOOOOE-Ol 

-1.900D00E+00 

-1.913723E+00 

-1.372324E-02 

0.72 

5.7000006-01 

-1 , 920  OQOE  +  OO 

-1.921479E+00 

-1.479117E-03 

0.08 

5.6000006-01 

-1.9500006+00 

-1.9436766+00 

6.324343E-03 

0.32 

15.9000006-01 

-2. 0000006+00 

-2.C03613E+00 

-3.6?  2798E-03 

0.18 

6.0000006-01 

-2. OOOOOOE+OO 

-1.9982456+00 

1.754805E-03 

0.09 

6.1000006-01 

-2. 0000006+00 

-1.994207E+00 

5 , 793437E-03 

0.29 

6.200000E-01 

-2.0000006+00 

-2. 0097776+00 

-9 . 777188E-03 

0.49 

6.3000006-01 

-2. OOOOOOE+OO 

-1.9924776+00 

7.5232596-03 

0.36 

6. 4cc.no  oE-oi 

-2. 000000E  +  00 

-2.0007196+00 

-7.189810E-04 

0.04 

6. 5000006-01 

-2.0000006+00 

-2.0061066+00 

-6.1078936-03 

0.31 

6.6000006-01 

-2.0000006+00 

-1.9913736+00 

8.6271736-03 

0.43 

6.7000006-01 

-2.000000E+00 

-2.005434E+00 

-5.4336456-03 

0.27 

o.dounoOE-Ol 

-2.0000006+00 

-2.001178E+00 

-1.177821E-03 

0.06 

o .  90  0  0006-01 

-2.(1000006+00 

-1.9931986+00 

6 .8020716-03 

0.34 

7. 000  0 006-01 

-2.000000E+00 

-2.007796E+00 

-7.795572E-03 

0.39 

Kill  scalf  error 

*  8.14 

•0. 70428555+00 
■0,58033315  +  00..... 

U.1D194725+00 ’  ’ 
0 . 936  J77t>g-Gl 
\o . 35898525-01 
•U-.  16460955-01 
*0.11679125-02 
u.70o0000E+0Q  _ 
«J.97b7<*326  +  00 
:0.l2t>2937E  +  00 
J.6C97020F+00 
•C. 24107255  +  00 
-0.D602277Ft.00 
■0.96045725+00  _ 
•0.76428555+00 
-0.48lD565Ft.0U 
•0.19139015-01 
0 . 536Q197F-U1 
0.8570564E-02 
0.67470065-02 
•0. 95705005-0.3"' 
0.70000005+00 
■0.32386415+00 
U.  15073715  +  00. 
u. 76741385-01 

■O.3466042E.-01 _ 

-0.1394463R-01 

0.6693835E-03__ 


■  6 . 26393925+01 
0.4223713E+00  • 

0 .1176537E+ 00  *« 
•0.7324159E-01  * 

■0.2126490E-01 
0.1611557E-01 
•0.1238799E-01 
0.1433546E+01_ 
0.7278940E+00  ■ 

0.1244368E+01. 
0.2066928E+00 
0.1139477E+00  ■ 

•0 . 8321966E+00  « 

,0.«  14272335+01 _ 

■0 . 26145716+01 
0 . 281 0  031E  +  0  0 
0.1114409E+00  * 

•0.15U3095E-01.  - 
•0.24178475-02  • 

■0.1191604E-01_  • 
■0.1226055E-01 
*0.3611200E  +  00._ 
0. 31533285  +  00  ■ 

0.3772702E-01.  * 
•0.7168264E-01  * 

l0,.2.1127.0.0.E^D.l _ 

0.1084986E-01 

_O..JLOOOOO.O£±.0.0_ 


C+d  our/fi+i') 


0.12665935+01 

0.74261345-01. 

0.1194020R+00 

0.1998052E-01 

0.37736895-01 

0.35945065-02 

0.13012545+0“.. 

0.1077931E+01 

0.4584138E+00 

0.10973115+01 

0.32308705+00 

0.7721945E+00 

0.12208595+01 

0.6542281E-01 

0 .52838505-01 

0.1778235E-01, 

0.11981195-01 

0.2507921E-02., 

0.33728155+00 

0.3513690E-01 

0.10707445+00 

C.9110935E-02 

0.3578409Er01_ 

0.2575038E-02 


0.1619213E-01 
-0.2387062E+0Q 
* 0.3952108E-01 
..  0.6323731E-01 
-0.1200203E-01 
-0.1260595E-01. 

_ 0..5816975E  +  00 

0.6039692E+00 

-C.1063192E+01. 

0 . 1551091E+00 
-0.6723846E+00 
-0.1173255E+01 


0 . 8896656E-02 
■0.1355646E+00 
•0.34537925-01 
0 . 97694  0  0E-02 
0.7475507E-02 
..0.H62288E-01 

.  0 . 1352903E-01 
•0.1964762E+00 
0.1921060E-01 
0 • 6247812E-01 
l0..  .93896555-02 
■0.4880430E-02 


1 


COURIER  S£ R I  PS  CURVE  FITTING 


:  =  5T  SITE  15/2  REVERSED  WITH  EXTENDED  FL-U._ _ 19  MARCH  1965 


A'i  0  1 

x 

.6435G0UE+01 

0.60000008*03 

A"PL.  cofff. 

ANGLE  (RADI 

N 

;35b2693F+0l 

-0.13859696*00 

1 

'.1097255E  +  01 

-C.1447923E+01 

2 

.  .1637494F+U1 

-0. 1122497E*00 

3 
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